matrix region=J(58, 8, 0)

*Ayacucho
scalar ir=16
matrix region[46,1]=1
matrix region[49,1]=1
while ir<=44{              /* begin loop */
matrix region[ir,1]=1
scalar ir=ir+1
}                         /* end loop */

*Centro
matrix region[15,2]=1
scalar ir=51
while ir<=57{              /* begin loop */
matrix region[ir,2]=1
scalar ir=ir+1
}                         /* end loop */

*Nor Oriente
matrix region[13,3]=1
matrix region[14,3]=1
matrix region[58,3]=1

*Sur Andino
matrix region[11,4]=1
matrix region[47,4]=1

*Huancavelica
matrix region[48,5]=1
matrix region[45,5]=1
matrix region[50,5]=1

*Lima-Callao
matrix region[2,6]=1
matrix region[8,6]=1
matrix region[9,6]=1

*Otros
* 1 3-7 10 12
scalar ir=1
while ir<=58{              /* begin loop */
matrix region[ir,7]=1-region[ir,1]-region[ir,2]-region[ir,3]-region[ir,4]-region[ir,5]-region[ir,6]
matrix region[ir,8]=1
scalar ir=ir+1
}                         /* end loop */




